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Abstract 

The two key phenomena occurring in the process of ceramic powder compaction are 
the progressive gain in cohesion and the increase of elastic stiffness, both related 
to the development of plastic deformation. The latter effect is an example of 
'elastoplastic coupling', in which the plastic flow affects the elastic properties of 
the material, and has been so far considered only within the framework of small 
strain assumption (mainly to describe elastic degradation in rock-like materials) , so 
that it remains completely unexplored for large strain. Therefore, a new finite strain 
generalization of elastoplastic coupling theory is given to describe the mechanical 
behaviour of materials evolving from a granular to a dense state. 

The correct account of elastoplastic coupling and of the specific characteristics 
of materials evolving from a loose to a dense state (for instance, nonlinear -or 
linear- dependence of the elastic part of the deformation on the forming pressure 
in the granular -or dense- state) makes the use of existing large strain formula- 
tions awkward, if even possible. Therfore, first, we have resorted to a very general 
setting allowing general transformations between work-conjugate stress and strain 
measures; second, we have introduced the multiplicative decomposition of the de- 
formation gradient and, third, employing isotropy and hyperelasticity of elastic 
response, we have obtained a relation between the Biot stress and its 'total' and 
'plastic' work-conjugate strain measure. This is a key result, since it allows an 
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immediate achievement of the rate elastoplastic constitutive equations. Knowing 
the general form of these equations, all the specific laws governing the behaviour 
of ceramic powders are finally introduced as generalizations of the small strain 
counterparts given in Part I of this paper. 

Keywords: Elastoplasticity; Large strains; Granular materials; Mechanical densifica- 
tion; Forming; Ceramic Materials. 

1 Introduction 

Mechanical cold compaction of ceramic powder involves the transition from a granular 
to a dense state. During this process, and strictly related to the development of perma- 
nent deformations, both the cohesion and elastic stiffness of the material increase. This 
occurs also under isostatic compression (which does not involve deviatoric strain) and 
is believed to be related at the microscale to the increase of the contact area between 
the grains (therefore, the effect should not be confused with a large strain effect). The 
increase of elastic stiffness with permanent deformation is a sort of 'inverse damage', 
which can be described making recourse to the concept of elastoplastic coupling, origi- 
narily invented to model elastic degradation, and employed in the Part I of this paper to 
describe the stiffening during plastic deformation of ceramic powder. However, a large 
strain formulation of elastoplastic coupling has never been attempted. That this formu- 
lation is not trivial can be deduced from the fact that elastic characteristics have been 
assumed to be independent of plastic deformation in all elastoplastic models proposed 
for soils (Borja and Tamagnini, 1998; Callari et al. 1998; Rouainia and Muir Wood, 
2000; Ortiz and Pandolfi, 2004) and in more general contexts (see among others: Simo 
and Miehe, 1992; Peric et al. 1992; Schieck and Stumpf, 1993; Simo and Meschke, 1993; 
Ibrahimbegovic, 1994). 

Since the existing large strain formulations do not appear to be easily generalizable to 
admit a coupling between elastic and plastic deformations, we have recurred to the early 
formulation by Hill and Rice (1973) (see also Hill, 1978; Petryk and Thermann 1985; 
Bigoni, 1996; 2000), which (although not explicitely mentioned) has been formulated in 
such a generality to include coupling. Following this approach, the level of generality is 
so high that the following choices are not required: stress/strain measures [except that 
these are work-conjugate (Hill, 1968)], elastic and plastic strain decomposition, elastic 
law, yield function, flow and hardening rules. After this framework is provided, the 
multiplicative strain decomposition of Lee (1969) and Willis (1969) is introduced. At this 
point, assuming that the elastic response be hyperelastic and isotropic we have proved 
that a general relation exists, in which the Biot stress is related to its work-conjugate 
'total' and 'plastic' strain measure (and to a generic set of scalar hidden variables). 
This achievement turns out to be crucial since it allows immediate use of the general 



formulation previously developed^. Finally, the coupling and hardening laws, the yield 
function and all constitutive relations provided in Part I of this paper for the infinitesimal 
theory are consistently generalized to include large strains. 

2 The skeleton of large strain elastoplasticity 
2.1 Some preliminaries on work conjugacy 

A broad constitutive framework for isothermal and time independent large elastoplastic 
deformations is presented, based on the concept of work coniugacy in the Hill sense 
(1968, 1978). In particular, employing Ogden's (1984) notation, a pair of symmetric, 
Lagrangean, stress and strain E^ measure^! are work-conjugate when the stress 
power density per unit volume in the reference configuration can be expressed as 

T (m) • E [m) = S-F, (1) 

where a dot over a symbol denotes material time derivative, F is the deformation gradient 
and S the first Piola-Kirchhoff stress tensor 

S = JTF~ T = KF T , (2) 

in which J = detF and T and K = JT are the Cauchy and Kirchhoff stresses, respec- 
tively. 

For integer (positive, null or negative) exponent m, we introduce the following La- 
grangean strain measures 

#M = —{U m - I), if m^O, 

m (3) 

E (0) = log U, if m = 0, 

where the logarithm of a tensor is defined as in (Ogden, 1984) and 

U = (F T F) 1/2 , (4) 

is the right stretch tensor. For a given m, is defined by ([3]) and the corresponding 
work-conjugate stress measure can be defined imposing eqn. (pQ). For instance, for 



1 The law between Biot stress and its work-conjugate strain measure could obviously be transformed 
into different (work-conjugate) stress/strain measures, but this would be cumbersome and useless, since 
the generality of the Hill and Rice (1973) approach allow us to use the obtained law directly. 

2 The notation T m = TT^T (or E m ) should not be confused with T (m) (or E (m) ). 

m times 



3 



m — 2, the Green-Lagrange strain results from eqn. ([3]) and the eqn. ([I]) provides for T^ 2 ' 
the second Piola-Kirchhoff stress tensors, 

E (2) = - (U 2 - I) , conjugate to T {2) = F X KF T . (5) 

A conjugate pair of stress and strain that will become useful later is formed by the Biot 
stress tensor and the strain measure E {1 \ defined as 

E {1) = U-I, conjugate to T (1) = l - (t {2) U + £/T (2) ) . (6) 

It is well-known however that it is not always the easy task of the two above examples 
to obtain the stress measure conjugated to a given strain of the form (j3]). For instance, 
the conjugate of the logarithmic strain E^ has a very complex form (Hoger, 1987), 
which simplifies to the rotated stress only when the two measures result coaxial, namely, 

E {0) = log U conjugate to T (0) = R T KR, (7) 

(where R = FU^ 1 is the rotation tensor) if and only if the following coaxiality condition 
holds true 

E (o) T (o) = T (o) E (o) ^ (logV)K = KhgV, (8) 

where V is the left stretch tensor, so that F = RU = VR (note that the above 
equivalence is an immediate consequence of the fact that the logarithmic function is 
isotropic). Condition OH]) is satisfied for isotropic elasticity, but may be not in more general 
contexts, such as for instance elastoplasticity (Sansour, 2001). It may be instructive for 
subsequent considerations to note that, when the coaxiality condition (JS} holds true, the 
following relation 

T (0) . E (°) = R . ( log yy (9) 

can be proved (Ogden, 1982) showing that the Eulerian stress and strain measures K 
and \ogV are work-conjugate (Hill, 1968). 



2.2 The basic assumptions of elastoplasticity 

Following Bigoni (2000), inelastic materials are considered that may at any stage of 
deformation exhibit a purely elastic response for appropriate loading. For these materials, 
elastic response is assumed to be a one-to-one relation between a work-conjugate pair 
T (m) and E {m) 

T (m) = f {m \E {m \}C), E {m) = £ (m) (T (m) ,/C), (10) 

where T and E are functionals depending on the prior history of inelastic deforma- 
tion through the unspecified set K, of variables of generic tensorial nature (thus embracing 
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scalars, vectors, second-order tensors and possibly higher-order tensors). For a purely 
elastic deformation rate, /C remains fixed, so that we have 

• (m) _ r • (m), • (m) _ „ r ' (m), , , 

T l = E[£T E K =M[T l % (11) 

where the fourth-order tensors E and M possess the minor symmetries induced by E and 
T, while the major symmetry is not a-priori requested (differently from Hill and Rice, 
1973). They are defined as 

f)f {m) f)F {m) 
E(£ (m) , JC) = — M(T (m) , JC) = ^=V^. (12) 
dE^ gj( m ) 

Tensors (IT2"|) obviously satisfy 

E = M^ X . (13) 
For an increment involving elastic and inelastic strain rates, we may write 

t (m) = E[E (m) ] - AE[P], E [m) = M[t (m) ] + AP, (14) 

where P G Sym sets the 'direction' (or the 'mode') of the irreversible deformation, which 
is given by 

„ i (m) „ * (in) 



The scalar A > appearing in eqn. (fl4j) is the plastic multiplier and vanishes for purely 
elastic response, namely, when JC = 0. 

A yield surface is assumed at each JC, which may be alternatively expressed in stress 
and strain spaces as 

f Tim) (T^ m \jC) < or as f E{m) (E^ m \jC) < 0, (16) 

thus defining regions of the T (m) or E {m) space, respectively, within which the response 
is elastic. 

2.3 Direct rate constitutive equations 

Prager's consistency condition requires 

f T (m) = f E (m) = 0, (17) 

when inelastic strain rate is different from zero. As a consequence, employing the stress 
space representation, the elastoplastic incremental constitutive equations can be written 

as 

E[E im) }--<Q.E[E im) }>E[P] if f T(m) (T^\jC) = 0, 
T (m) = { 9 ( 18 ) 

E[E im) ] if / T(m) (T (m) ,/C) < 0, 



where the operator < ■ > denotes the Macaulay brackets, i.e. V a £ R, < a >= 
(a + |a|)/2. Moreover, the symmetric second-order tensor 



Q 




(19) 



is the yield function gradient and the plastic modulus 



g = h + Q-E[P] 



(20) 



is assumed to be strictly positive (a negative plastic modulus would correspond to a so- 
called locking material, not considered here). In the Hill (1967) notation, the hardening 
modulus h in (1201) describes: 

• hardening when positive, 

• softening when negative, 

• perfect plasticity when null. 
The hardening modulus is defined as 



and, as Hill (1967) remarks, hardening and softening are not measure-invariant concepts, 
in the sense that h depends on the choice of T*™- 1 and E^ m \ Therefore, the above 
nomenclature is, to some extent, arbitrary. Moreover, we remark that, in addition to 
h, also Q, P and E are measure-dependent. On the contrary, the plastic modulus g 
can be shown to be measure- independent (Hill, 1967; Petryk, 2000). Note also that 
all quantities appearing in the rate equations ( !T8|) fully depend on the entire path of 
deformation reckoned from some ground state. 

2.4 Inverse rate constitutive equations 

Under the assumption of positive hardening, the rate equations (|T8|) can be inverted 
to relate the material derivative of the strain measure to the material derivative of the 
work-conjugate stress. In particular, taking the scalar product of the first equation in 
(USD with Q gives 



Ah 




•K, 



(21) 



dlC 



QT 



Q • E[E Kni '] 



•(m), Q-E[P] 



< Q ■ E[E (m> ] > . 



(22) 
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In the case when h > 0, we note that 



sgn(Q.E[ J E; (mj ]) = sgn(Q.t imj ). 
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Therefore, assuming positive hardening, h > 0, and using fl22|) . we obtain the inverse 
constitutive equations 



The rate constitutive equations (TTSJ) or (I2"3"j) represent a broad constitutive framework, 
within which 

all possible choices of and are equivalent and the requirement of 
material frame indifference (Truesdell and Noll, 1965) is never violated. 

Moreover, the framework is so general that it does not imply any particular choice of 

• elastic and plastic strain decomposition, 

• hypo- or hyper- elastic law, 

• yield function, flow and hardening rules. 

It is however clear that in order to set up the constitutive modelling of a particular 
material, we need to introduce specific laws. This objective will be pursued in three 
steps of deceasing generality in the following: first, we will introduce the multiplicative 
elastic and plastic strain decomposition and requirement of objectivity and isotropy of the 
elastic constitutive law; second, a form of elastic constitutive equation will be proposed, 
depending on plastic deformation and thus capable of describing the elastic behaviour 
of granular and dense materials; third, yield function and hardening laws are introduced 
as simple generalizations of the rules formulated in Part I of this paper under the small 
strain assumption. 

3 The multiplicative decomposition and the elastic 
law 

The multiplicative decomposition of deformation gradient F into elastic F e and plastic 
F p components introduced by Lee (1969) and Willis (1969) is adopted (Fig. [Q. 



According to eqn. (T21j) . using the left polar decomposition F = VR, we introduce the 
elastic and plastic left stretch and rotation tensors V e , V p , R e and R p , satisfying 




(23) 



F = F e F p . 



(24) 



F 

- 1 - e 



V R F 



V R 



(25) 
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deformed 




configuration 



Figure 1: Reference, current and intermediate configurations. 



while using the right polar decomposition F = RU, the right elastic and plastic stretch 
tensors U e and U p are defined so that they satisfy 

F e = R e U e , F p = RpU p . (26) 

A crucial expedient, employed also by Ortiz and Pandolfi (2004), to describe the 
behaviour of granular materials is to refer to the logarithmic strains defined as 

e = log V, e e = log V e , e p = log V p , (27) 

and 

£ {0) = log U, = log U e , Ef = log Up. (28) 

The interest in employing definitions ( 1271) and ( 1281) is that these allow a decoupling 
between the volumetric logarithmic elastic and plastic deformations, namely 

tr e = tr e e + tr e p = tr E {0) = tr E^ + tr Ef , (29) 

which, employing the usual definition of J and noting the property (for every symmetric 
tensor A) 

tr(log A) = log(det A) , (30) 

can be written as 

log J = log J e + log J p. (31) 



4 Objectivity and isotropy of elastic response 



We refer now to an isotropic elastic law relating the Kirchhoff stress K to the elastic 
deformation gradient F e in the generic form 

K = K(F e ,k i ), (32) 

where function K may depend also on generic plastic scalar variables ki, assumed in- 
variant with respect to every symmetry group of the material and change in observer. 

In general, the elastic response must be objective, but in addition we assume for 
simplicity that the elastic response be isotropic. The latter requirement implies the 
coaxiality condition (jSJ), a requisite more important than it may appear, since it ensures 
the work coniugacy (|7j) and its Eulerian counterpart. Therefore, function K is assumed 
to satisfy: 

PI. The objectivity requirement 

K(F e , h) = R T K(RF e , ki)R, \/R G Orth+. (33) 

P2. The isotropy requirement 

K(F e , ki) = K(F e R, hi), \/R G Orth+. (34) 

Though the generalization of the formulation to anisotropy of the elastic response is defi- 
nitely important to capture certain experimental evidence connected to the development 
of various form of instabilities [see Gajo et al. (2004) for a discussion relative to the 
infinitesimal theory in the context of granular media], we note that this generalization 
is for the moment lacking for granular material subject to large strains, even in the 
relatively simple setting in which cohesion and coupling are neglected. 

As a consequence of isotropy, eqn. fl34l) . the rotation in the left polar decomposition 
does not alter the values of function K, 

K(F e ,ki) = k(V e ,ki), (35) 

so that function k depends only on the elastic left stretch tensor. Noting the identity 

F e = RUU; l R T p , (36) 

isotropy and objectivity allow us to introduce the following transformations 

K = K(RUU p l R T p , ki) = k{RUU~\ h) = RK{UU~\ h)R 7 \ (37) 

so that we get 

K = K(F e , h) = RK{UU-\ ki)R T . (38) 
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Employing now the rotated stress R T KR, the constitutive law (132]) can be written in 
the form 

R T KR = k(UU~\ hi), (39) 

relating the rotated stress to the total and plastic right stretch tensors. 

Now, the rotated stress is related to the Biot stress through (Ogden, 1984) 

T {1) = - (U~ 1 R T KR + R T KRU l ) , (40) 
so that in conclusion we obtain the elastic constitutive law in the form 

r« = i (u- x k{uu;\ h) + k(uu~\ k^u- 1 ) , (4i) 

relating the Biot stress to the global and plastic right stretch tensors. 

Since eqn. ([6]) shows that tensor E^ 1 ' is the right stretch tensor U with the identity 
subtracted, 

eqn. (T4TT) expresses a relation between the two work-conjugate measures 
and E^ of the type (fit)]) , in which the set /C is now including U~ l and the 
scalar s k{. 



5 Formulation of the rate model 

Until this point, all the equations are at a high level of generality; now, to develop the 
model for powder densification, further specific laws are introduced, including the par- 
ticular hyperelastic-plastic coupling rule, yield function and hardening laws. Since these 
laws are essentially extensions of those already employed in the small strain formula- 
tion, details on physical motivations determining the specific choices will be omitted for 
conciseness. 

5.1 Elastoplastic coupling 

The elastic properties of granular materials can be described by a hyperelastic nonlinear 
law providing a generalization to finite strains of the corresponding equation introduced 
in Part I, Section 2.5. This generalization is represented by the following potential 

0(e e ,e p ) = -^(tre e ) 2 + ctre e (42) 
+ (po + c) 

where k is the elastic logarithmic bulk modulus, po is the initial confining pressure, c, d, 
and are scalar parameters depending on the volumetric plastic strain tr e p = tr Ey , 



t r e e 

/«K, 



+ fjL€ e 
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providing the elastoplastic coupling. The dependence of parameters c, d, and /i on the 
volumetric plastic strain is made explicit by the following equations 



li + c ( d- - ) //i, 



d = 1 + B < p c -p cb >, 



c = Coo [1 - exp (-r < Vc ~ Pcb >)] , 



and 



exp (tiE^^) — 1 = — ai exp f -J — a 2 exp ( -J 

V p / V Vc) \ Vc) 



(43) 



(44) 



where /x , /ii, -B, p c fe, Oxd, T, a 1( a 2 , A 1; and A 2 are positive material constants. It can be 
noted that eqns. (143]) do not include 'geometrical terms' and thus coincide with those of 
the small strain formulation, whereas eqn. (14"4"|) has been consistently generalized. 

The potential fT421) represents a isotropic function of the logarithmic strain and can be 
written as the sum of a volumetric and a deviatoric component. The deviatoric potential 
coincides with that employed by Ortiz and Pandolfi (2004) in the special case of null 
cohesion c = and null coupling d — 1. 

The Kirchhoff stress can be obtained from the potential f T42|) as 



so that it results in the form 



K 



■-/itre e + c+ (p + c) 



K 



<9e e ' 
tr e f 



exp 



tr e f 



(45) 



J + 2/xe e . (46) 



Eqn. ( l4"6"j) implies that the Kirchhoff stress and the logarithmic strain are coaxial, so that 
these become in the present context work-conjugate stress and strain measures, eqn. 0. 
The elastic constitutive law (H6l) can be written in the form fl32l) with 



K(F e ,ki) 



--/xtrlog F e F 2 e +c 



(47) 



+ {Vo + c) 



d _l\ tr log F e F J e 



d 



2k 



exp 



trlogF e F^ 



2d x ' n 



K 



+/ilog F e Fl. 

where the set ki is now including c, d, and /x. Obviously, eqn. (14T|) can be expressed in 
the form (jHJ), not reported for conciseness. 

We are now in a position to write down the fourth-order elastic tensor E defined by 
eqn. ( I12[) i. This takes the form 



E = -- ( U- 1 ® KU 1 + KC/" 1 ® C/" 1 



(4J 



+ 



c7X 



(Jit/; 1 ) 



x=uu~ 
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so that, employing eqn. f|471) . we obtain 



dK(X,ki) 



dX 



(49) 



r — T 1 f dlogY 



Y=XX' 



where 



and 



K,(X) 



Po + c 
2k 



d--, + d~ 1/n 
a 



exp 



trlogXX J 

2dV n R 



d log Y 
dY 



00 I 1 \n+l n ~ 1 



Y-I) 



(y - iy 



n—l—r 



(50) 



(51) 



n=l r=0 

a formula that can be found in (Hoger, 1987) or deduced from (Itskov and Aksel, 2002). 
Note that four tensorial products between second-order tensors A and B have been 
employed, which can be defined, with reference to every tensor C, as 

(A ® B) [C] = (C • B T )A, ( A ® S) [C] = I A (C + C T ) B T , 

2 (52) 

(A®B)[C} = ACB T , 
so that the following property holds 



(A®B)[C] = AC £? 



(53) 



5.2 The irreversible strain rate 



Eqn. (U5l) i defines the irreversible strain rate, which, accordingly, may be calculated 
taking the derivatives of eqn. ( 14 ip in which K is given by eqn. (T47j). We obtain 



AP = G 



(0)1 



(54) 
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where tensor G, assumed positive definite, is given by 



G = — ^E- 1 
2 



oc oc 



O/I Ojl 



+ 2 E OX + dX~ 



x=uu- x 



1N dexpM 0) 



in which 



<9K 

"ft 7 



1\ trlogC/t/; 2 f/ 

~ d 2k 



— cxp 



trlogC/C/r 2 C7 



2d l ' n 



K 



dK _ \(po + c) tr log UU~ 2 U 
~dd ~ 



2k 



1 + T> 



1 



trloe UU~ 2 U 



2nrf 1+1 /' 



exp 



2d l ' n 



K 



and 



dK 



-- tr log UTJ-/U / + log UU; Z U, 



6 = 



<9/z V 3 



Coo r g(p e - p cfc ) exp [-Tjpc - Pcb)} P 2 C exp(tr E^) 

7 m " r - ^ 

aiAi exp + a 2 A 2 exp 

. £#(p c -p cfc )p^xp(trgf ) 

3 " 4 f AA . / A 2 
aiAi exp + a 2 A 2 exp 

\ PcJ \ p c 
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U = [d-\)^ 2 + c(l + ^)^ 3 . (61) 



d J ™ V d 2 

Note that the exponential of a tensor [see e.g. (Itskov and Aksel, 2002) for the definition] 
has been introduced in eqn. (|55|) . together with its gradient, defined as 

UJ ^p „=l r =o 

a formula given by Itskov and Aksel (2002). 
5.3 The yield function 

In the absence of ad hoc experimental results, we employ for simplicity a yield function 
with same form adopted for the infinitesimal theory (see Part I of this paper and Bigoni 
and Piccolroaz, 2004), where now the Cauchy stress is replaced by the Biot stress 
This can be pursued re-defining the invariants p, q and 9 in terms of Biot stress 

trTW r— . 1 _! (3y/3 J 3 \ , R „, 

P= 3 — , g=V3^2, 9 = -cos I 2 ^ 3/2 1 , (63) 

where 9 G [0,7r/3] and 

J 2 = -devT« -devT^ J 3 = -tr (devT^) 3 , 

2 trT« 3 (64) 

3 

As a consequence, the yield function takes a form of the type f T m (T^\ fC) < 0, namely, 

F(T^,p c ,c) = f(p,p c ,c) + ^-, (65) 

where p c and c are the parameters governing the change in shape of the yield surface 
caused by the hardening (as in the infinitesimal theory), and 



-Mp cV /($ - $ m ) [2(1 - a)$ + a] if $ G [0, 1], 
f(p,p c ,c) = { (66) 
+oo if $£[0,1], 



in which 

p + c(tr M 0) 



$ = — (67) 

p c {tiE^) + c{tvEf) 
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and 



9(0) 



cos 



1 



7T 1 

~ o cos_1 (7cos36») 

6 



(68) 



Note that M, m, a, (3, and 7 are material parameters with the same meaning as in 
the infinitesimal theory (already described in Part. I of this paper and by Bigoni and 
Piccolroaz, 2004). 

The yield function gradient 

dF(T^,p c ,c) 



Q 



dT w 



(69) 



can be obtained directly from the yield function (|66|) or from the corresponding equations 
in Part I of this paper (Appendix A), with replacing cr. 



5.4 Flow rule and hardening modulus 

The flow mode tensor P is postulated in the form 

P = Q-^e(l-*)I, (70) 

where < e < 1 is a nonassociativity parameter, null for associative flow rule. 
The hardening modulus can be calculated from the definition (1211) in the form 

• (dF. dF \ 

Ah = -[ — V c + -5-c , (71) 



where 



\ dp c dc 



— = -iWV (<£> - $ m ) [2(1 - a)$ + a] (72) 
op c 

+M Pc(P + c ) ( X - m^ m - 1 ) [2(1 - a)$ + a] + 2(1 - a) ($ - $ m ) 



(p c + c) 2 2v/($ - $ m ) [2(l-a)$ + 



and 



dF _ PciPc-p) (73) 
(1 - m^™- 1 ) [2(1 - a)$ + a] + 2(1 - a) ($ - $ m ) 



2^($ - $m) [2(1 - a )$ + 



in which 



pgexp(tr4°)) (o) 

Pc = "t: — — ; — p ' ( } 

aiAi exp + a 2 A 2 exp 



Vc J \ Vc 
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and 

c = c^r H(p c - p cb ) exp [-r(p c - p ch )) p c . (75) 

Parameters A 1; A 2 , 5i, a 2 , c^, T, and p c b have been introduced and motivated in Part I 
of this paper. 



5.5 The rate constitutive equations for the compaction model 



The elastoplastic incremental constitutive equations ( fl8l) . written in terms of Biot stress 
T^ 1 ' and conjugate strain E^ 1 ' = U — I, take the form 

[ E[E (1) ] - - < Q • E[E {1) ] > E[P] if F(T« Pc ,c) = 0, 
t (1) = I 9 ( 76 ) 

i E[£ (1) ] if F(T«,p c ,c)<0, 

where the elastic tensor E is given by eqns. fj48|) - fl5T|) . the yield function F(T^\p c ,c) 
by eqns. fl65|) -f l68|) and the yield function gradient Q and flow mode tensor P by eqns. 
fl69|) and f JTOj) . The plastic modulus g is provided by eqn. f l20|) . where /i, the hardening 
modulus, is obtained substituting eqn. f|51|) into eqns. f !71|) - fl75]) and eqn. fITT]) . thus 
yielding 

h = -( t^Pc + -t^c , (77) 



where dF/dp c and dF/dc are specified by eqns. (|72]) - (175|) and 

p2 ex p(trM 0) ) 1r , 

«= = -- — 7 AA 7 M trG [F1 - 

aiAi exp I I + a 2 A 2 exp I -— I ( 7 j 

c = c^r ii(p c - p cb ) exp [-r(p c - p c6 )] p c , 

in which tensor G is explicited by eqns. (!55|) -( l62|) . 



6 Conclusions 

A new, consistent generalization to large strains of elastoplasticity theory with coupling 
between elastic and plastic properties has been given, based on work-conjugate variables 
and isotropy of the elastic response. This has permitted the extension to large strains 
of the model introduced in Part I of this paper, to describe granular materials becoming 
cohesive during mechanical, cold densification. Therefore, the model allows the simula- 
tion of forming processes of green bodies from ceramic powders, including situations in 
which large strains are involved. 
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